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ABSTRACT 

The source of hot gas in elliptical galaxies is thought to be due to stellar mass 
loss, with contributions from supernova events and possibly from infall from a 
surrounding environment. This picture predicts supersolar values for the metal- 
licity of the gas toward the inner part of the galaxy, which can be tested by 
measuring the gas phase abundances. We use high-quality data for 10 nearby 
early-type galaxy from XMM-Newton, featuring both the EPIC and the Reflec- 
tion Grating Spectrometer, where the strongest emission lines are detected with 
little blending; some Chandra data are also used. We find excellent consistency 
in the elemental abundances between the different XMM instruments and good 
consistency with Chandra. Differences in abundances with aperture size and 
model complexity are examined, but large differences rarely occur. For a two- 
temperature thermal model plus a point source contribution, the median Fe and 
O abundances are 0.86 and 0.44 of the Solar value, while Si and Mg abundances 
are similar to that for Fe. This is similar to stellar abundances for these galaxies 
but supernovae were expected to enhance the gas phase abundances considerably, 
which is not observed. 

Subject headings: cooling flows - galaxies: abundances - galaxies: elliptical and 
lenticular, cD - galaxies: individual (NGC 720, NGC 1399, NGC 3923, NGC 
4406, NGC 4472, NGC 4553, NGC 4636, NGC 4649, NGC 5044, IC 1459) - 
X-rays: galaxies 



1. Introduction 



Early-type galaxies possess an interstellar medium that is dominated by hot gas (3-10 x 
10 6 K), although the mass of gas can vary widely between systems. The origin of the hot 
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gas is not entirely a settled issue, but it is probably the result of mass loss from stars within 
the galaxy as well as infall onto the galaxy, especially when it lies in a galaxy group. The 
abundances of this gas reflect its history and can potentially inform us as to the number of 
Type I and Type II supernova that must have been present. There have been a number of 
surprises in the abundance measurements, such as that the abundance is lower than initiall y 



predicted for gas shed from stars and enhanced by supernovae (e.g., lArimoto et al.l 119971 ). 
Also, the values for the gas abundances have varied considerably, for different models applied 
to the same galaxy, and between galaxies, so a uniform picture has been slow to emerge. 

There are a variety of issues that face investigators when determining abundances within 
early-type galaxies (or other systems with thermal gas). There is the problem of instrumental 
calibration, which can be notoriously difficult, despite dedicated efforts by the scientific staff. 
It is often several years after launch of a mission before most of the important calibration 
issues are understood. As the calibrations for XMM and Chandra have matured, this is a 
good time to examine the spectra of similar objects and compare the results. 

Another issue is that the derived metallicities a re sensitive to the num ber of spectral 



components used in a model. This was pointed out by lTrinchieri et al.l (119941 ). among others, 
who showed that when a single-temperature thermal model was applied to the luminous 
emission from an elliptical galaxy, the derived metallicity was significantly lower than when 
a two-temperature model was used. Often, the \ 2 is acceptable for both models, so without 
further information, it is difficult to identify the correct model. The problem is one of 
resolution, both spatially and spectrally. On one hand, the temperature can vary within 
a galaxy, so by analyzing projected emission, there are multiple temperature components. 
Also, there are point sources, mainly due to low-mass X-ray binaries, and these provide a hard 
continuum that must be accounted for. Even with Chandra, not all of the individual point 
sources ca n be excluded, alt hough their collective spectra are fairly constant from galaxy 



to galaxy (llrwin et al.l 120031 ). so modeling of this component is tractable. Spectrally, one 
could identify the need for various spectral components if it were possible to measure lines 
of various ionization states for the same element. ROSAT did not have sufficient spectral 
resolution nor bandpass coverage to constrain many of the important parameters nor did it 
have the spatial resolution to remove point sources. The ASCA satellite could measure the 
high energy contribution from the X-ray binaries, but had very poor spatial resolution and 
there were calibration issues at the important low-energy part of the detector. 

Some of these issues are resolved by using Chandra and XMM observations. The Chan- 
dra data have excellent spatial resolution, so most point sources can be excluded, and the 
combination of spectral resolution and calibration is superior to its predecessors. In com- 
parison, XMM has poorer spatial resolution but more collecting area and a relatively high 
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dispersion grating spectrum is obtained for all on-axis targets. 



There is a range of results that seems puzzling, as both high and low metallicities are 
found in optically similar galaxies. Individual XM M observations show a sim ilar range of 
behav ior, with subsolar ab undances (referenced to lAnders fe Grevessd Il989l). as in NGC 
6251 JSambruna et al.ll2004h . N GC 3585, 4494 and 5322 Jo'Sullivan fe Ponmanl l2004f). near- 
solar value s, as in NGC 4649 (jRandall et al.ll2006l ). through the supersolar values seen in 
NGC 507 (IKim fe Fabbianol I2004J ) . Similarly, indivi dual Chandra obervatio ns of galaxies 
can be of near solar metallicity, such as in NGC 1316 (IKim fe Fabbianoll2003l ) or NGC 4649 
( Randall et al.ll2004f) . but other galaxies can show quite low abundances, such as NGC 1291 
( llrwin et al.ll2002l ). Two surveys of galaxies with Chandra also show a range metallicities, 
but with no meaningful correlation between the stellar metalli cities and the metallicities of 
the X-ray emitting gas (lAtheyl 120071 ; Humphrey fe Buotell2006l ). 



If there is a consistent trend, it is with the ratios of some of the elements, such as in the 
XMM observation of NGC 4649, where the light elements (N, C, O, Ca; this is dominated 
by O) are 0.22 of the Solar value, Fe is 0.7 Sol ar, while the combi nation of Ne, Mg, Si, S, Ar, 
Al, Na, and Ni is 1.7 times the Solar value (jRandall et al.l 120061 ) . These results of O being 
low relative to Fe are very common, as are the relatively hig her values of Mg and Si. S uch 
ratios may be in conflict with expectations from supernovae ([Humphrey fe Buotdl2006l ). 



The above studies either have been in-depth investigations of a single galaxy with a 
single instrument (often with different models and approaches), or a study of several galaxies 
with a single instrument, but there has not been a comparison of results across instruments. 
One of our objectives is to analyze several galaxies using both Chandra and XMM data 
and applying the same models to both data sets. This helps to address whether there are 
serious calibration differences between instruments that might yield misleading results. Also, 
there has been little use of the Reflection Grating Spectrograph (RGS) data from XMM, 
which has the highest spectral resolution and where individual strong spectral lines can be 
detected and identified. With advances in how to model extended sources observed with the 
RGS, these spectra become extremely useful and we present these for several bright sources. 
Finally, we examine the differences in the abundances that result when different instrumental 
backgrounds are adopted. 
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2. Sample Selection and Data Reduction 
2.1. Sample Selection 



The sample is based on the work of iBrown fc Bregmanl (119981 ). who obtained a complete 



sample of 34 early-type galaxies which were optically selected and flux-limited. We selected 
10 sources in their sample which were observed by XMM-Newton, and are bright enough to 
use both of the data of the European Photon Imaging Camera (EPIC) and the Reflection 
Grating Spectrometer (RGS). This is mainly determined by the quality of the RGS spectra, 
which were included in our sample by presenting enough emission line features in the preview 
spectra. Several of these galaxies also have high-quality Chandra data, which were used to 
determine abundances as well. The observations that were used in this paper are given in 
Table [U 



2.2. Chandra ACIS Data Processing 

Chandra observations of NGC 4649, NGC 4472, and NGC 1399 were obtained from the 
HEASARC online archive ( |http : / /hea sar c . gsf c . nasa. gov /docs/ archive . html | . The data for 
each observation were processed in a uniform manner following the Chandra data reduction 
threads using ciao v3.4. Times of high background were removed from the data, and all 
images created were corrected for exposure and vignetting. Point sources were identified 
and subsequently removed using "Mexican-Hat" wavelet detection routine wavdetect in 
CIAO on the 0.3-6.0 keV band image. For each galaxy, we extracted spectra from a region 
1' in diameter centered on the nucleus of the galaxy (the same region as was used with 
XMM-Newton data for these galaxies). Background spectra were extracted from a region 
near the edge of the S3 chip for each galaxy. While the chosen background region may 
contain extended emission from the galaxy (which typically fills the entire S3 chip), it only 
represents ~1% of the total emission in the 1' diameter aperture, so precise treatment of the 
background is unnecessary. The spectra were grouped such that each spectral bin contained 
at least 25 counts. 



2.3. XMM-Newton EPIC data reduction 

The EPIC and RGS data of the 10 sources were processed with standard procedures 
using XMM-Newton Science Analysis System (SAS, Version 7.0.0). 

For the EPIC data, high background periods (counts rate >0.35 cnt s _1 for MOS CCDs 
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and >1 cnt s _1 for PN CCD) were removed by examining light curves of events with photon 
energy greater than 10 keV and PATTERN = 0, which is very sensitive to the X-ray 
flares. Photon events with PATTERN < 12 were retained for the MOS cameras while 
PATTERN < 4 were retained for the PN camera. We chose the photon events in the 
energy range 0.3-7.0 keV for the MOS and 0.4-7.0 keV for the PN to eliminate the calibration 
uncertainty at the lower energy band (<0.3 keV for MOS and <0.4 keV for PN) and the 
uncertainty at the high energy band (>7.0 keV) due to very low photon counts. 

Two different source regions were extracted from the cleaned event lists. The regions 
are 1' and 4' in diameter, centered at the peak of the source X-ray emission, which is also 
the optical center of the galaxy. The reason for this is that the smaller region will be more 
uniform than the larger region, when radial gradients exist in the gas, such as temperture 
gradient. The larger region is representative of the average properties of the entire galaxy, 
so this analysis might be useful for comparisons with more distant galaxies. The larger 
diameter region will sample a broader variation of gas properties, such as in temperature 
and abundance, so by comparing the spectral fitting results from the two regions, we can 
assess the implied variations. 

We subtract the background regions on the same observation fields, which are a few 
arc-minutes away from the central sources and contain little hot gas emission from the 
central sources; point sources were excluded. The background regions are larger than the 
source regions to compensate the vignetting effect on the edge, and are circular regions with 
diameters of 1.67' and 5.83' for source region of 1' and 4', respectively. Redistribution matrix 
files (RMF), which contain the instrumental responses, and ancillary response files (ARF), 
which contain the effective area of the detector, were generated with SAS. The spectra were 
grouped so that each energy bin contains a minimum of 25 counts, appropriate for \ 2 fitting. 



2.4. XMM-Newton RGS data reduction 

For the RGS data, high background periods were screened by examining the light curve 
of CCD9 (excluding the source region), which is the most sensitive CCD to be affected by 



the b ackground flares. The screening criteria is the same as employed by iTamura et al. 



( 120031 ) . because we are using the same template RGS background. Events were selected in 
the energy range of 0.33 - 2.48 keV (5 - 38 A), and the first order spectra of RGS1 and 
RGS2 were extracted. 



The RGS is a slitless spectrograph with a field of view of 5' in the cross-dispersion 
direction and 1° in dispersion direction (covering the whole diameter of MOS field of view). 
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The RGS1 and RGS2 are perpendicular to each other and therefore sample different projected 
spatial regions of the galaxies. The size of the extended emission can be inferred from the 
their EPIC images, and we can fit the suface brightness distribution with a "beta" model, 
S x oc (1 + (r/r c ) 2 )- 3 ^ +1 / 2 , where (3 ~ 0.5 and r c < 0.2'. Although some of these galaxies 
are detectable in X-rays to 5' the flux is concentrated. Based on the surface brightness 
distribution, we defined two source extraction regions in the cross- dispersion direction. For 
a cross-dispersion width of V (to be compared with 1' diameter regions), 90% of the flux is 
extracted, while in our larger extraction width of 3.5', 99% of the flux is extracted. These 
extraction regions are somewhat different than the circular apertures used to analyze the 
EPIC data, which is a consideration in the spectral fitting process. If there is a gradient 
in spectral properties, the narrower extraction width will suffer less spectral blending. The 
advantage of the spectrum with the larger extraction width is higher S/N, provided that the 
source remains brighter than the background. 

When extracting a spectrum, a background must be removed and the standard pro- 
cedure is to use a local background, defined as a region in the same CCD, after exclud- 
ing sources. However, some of these galaxies fill the field of view of the RGS CCDs, so 
there is no emission-free background region that can be used. Therefore, we used a back- 
ground from other ob servations. For our observations, we used the RGS template background 



( ITamura et al.ll2003l ). combined with several other blank sky observations. The background 
is selected with the same region on RGS CCDs of the template as the source region. Another 
advantage of this method is that the background is better defined than a local background 
due to the very long exposure times of the background fields. 

Following screening and background subtraction, the RGS spectra were binned so that 
each energy bin contain a minimum of 25 counts, making the data appropriate for \ 2 fitting. 



3. Spectral Fitting Techniques and Results 

3.1. Spectral Fitting Techniques 

Originally, we anticipated using the RGS spectral fits to yield unambiguous abundance 
results, as individual lines are resolved from OVIII (19 A) and highly ionized Fe, such as 
FeXVII (15, 17 A) and FeXVIII (16 A) (see Fig. Q]), whereas they are badly blended in 
the EPIC data. These RGS data provide clear identification of the line strengths and the 
absence of certain important lines, such as the OVII triplet (20.6-21.1 A). However, the 
absolute abundance is a ratio of the line fluxes to the continuum fluxes, and this is less well 
determined in the RGS data (compared to the EPIC data) for two reasons. The first is that 
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the S/N is poorer in the RGS data due to fewer photons. The second point is that there is a 
power-law continuum contribution from LMXBs, which can be constrained at high energies 
from the higher S/N EPIC data. This power-law component cannot be properly constrained 
with the RGS data because of lower S/N and because the instrument has no response above 
2.5 keV, where the power-law is most distinct. Therefore, we found it useful to consider the 
EPIC and RGS spectral fits separately as well as jointly. All spectral fitting was performed 
within XSPEC (Version 11.3.2). 

As discussed below, there is good consistency when fitting the EPIC instruments sep- 
arately, so we permfom a simultaneous fit of a single model to the MOS1, MOS2, and the 
PN data with their normalizations free to fit. For the the rmal emission from thermal gas, 
we used the Astrophysical Plasma Emission code (APEC, ISmith et al.ll200ll ). which mod- 
els optically-thin emission from an isothermal plasma of some abundanc es (vapec) . The 
optically-thin assumption should be appropriate in most cases, although IXu et al.l (120021 ) 
demonstrate that optical depth effects occur in the strongest Fe XVII lines for the galaxy 
NGC 4636. Variations in temperature are accomodated by introducing two vapec models 
with different temperatures but same set of abundances. For the fit s, we fix the X-ray 
absorbing column at the Galactic HI value (jDickey fe Lockmanl Il990l ). and we fit for the 
temperture, the abundances of Fe, O, Si, and Mg, and the normalization. When the data 
had sufficiently high S/N, we also fit for Ne, Ni, and S, otherwise, they were tied to Fe with 
the ratios given by the Solar values. All other elemental abundances were tied to Fe, except 
for He, which was set at the Solar value. The abundances are referenced to the Solar values, 



which are defined he re by the values of iGrevesse fe Sauvall (119981 ). In many earlier works, the 
Solar abundances of Anders fe Grevessel ( 19891 ) were used, which have higher concentrations 
of O and Fe by 26% and 48%. Using the newer Solar abunandances, the relative abunances 
are raised. 

When fitting for the power-law component, most point sources are unresolved, unlike 
the situation for the Chandra observations. Based on Chandra observations of point sources 
in early-type galaxies, the summed x-ray point source s pectra can be rep resented by a power 
law model with an average photon index T of 1.56 (jlrwin et al.ll2002l ). This component 
dominates mainly the hard x-ray band (> 2 keV), but also contributes some emission to the 
soft x-ray band (< 2 keV). In the models, the contribution from this model is implemented 
by including a power-law component, with a fixed photon index of T = 1.56 and a free 
normalization; this component adds one free parameter. From our tests, we find nearly the 
same value for the power-law normalization when we only consider data above 3 keV (and 
no thermal component) or when we fit the entire energy range of the spectrum, using a 
power-law and a vapec component. 
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The procedure for fitting the Chandra spectra are similar except that most of the point 
sources are excluded from the spectral fits. Still, there remains the contribution from un- 
resolved point sources, so we included a power-law component, with fixed T to 1.56, along 
with the themal component (vapec) and a Galactic absorption term (wabs). 

For the fits to the RGS data, there are a few other considerations, the most important 
being the broadening of the lines along the dispersion direction due to the extended nature 
of the emission. A correction for this can be made if the spatial distribution of the emission 
is known, and for this purpose, we used the EPIC MOS1 image (equivalently, the other 
images could have been used). Using this image, the broadening correction is implemented 
by the task rgsxsrc within XSPEC. In the first set of fits to the RGS data, we fix the ratio of 
the luminosities from the thermal plasma component [vapec) to the power-law component 
(powerlaw) at the value given by the fits to the EPIC data. For the thermal plasma model, 
we use the RGS data to fit for the temperature, the normalization, and the abundances for 
O and Fe, which are the element abundances that are best constrained by the RGS data. 
Abundances for Mg and Si are determined whenever possible, otherwise they are tied to Fe. 
In addition, lines are detected from Ne and N for some sources, but the S/N is not sufficient 
for useful determinations, so the abundances of N, Ne, S, and Ni are tied to Fe. 



3.2. Comparison of ACIS-S and EPIC Spectral Fits for Bright Sources 

In the comparison between the four imaging devices (see Fig. [2]), the results for the inner 
1' diameter of NGC 4649 were fairly consistent, especially between the different EPIC instru- 
ments. Using the model described above, best fits were obtained separately for the ACIS-S3 
data, the PN data, the combined MOS1+MOS2 data, and the combined RGS1+RGS2 data. 
This galaxy is an important case because the temperature gradient is small and there is 
no evidence for a multi-temperature gaseous medium. We find that the temperatures were 
nearly identical for all detectors, but there were differences in the abundances of O and S 
between the XMM EPIC instruments and the Chandra ACIS-S detector. The S abundance 
obtained from the ACIS data is 75% lower than that derived from the EPIC data, and the 
90% error bars do not overlap. However, the two abundance measurement differ at the 1.8a 
level, so this is not a major discrepancy. A similar, but slightly greater discrepancy exist 
for the O abundance, where the EPIC value is about 2a above the ACIS value. Aside from 
these differences, the abundances follow rather closely for the different instruments in that, 
relative to Solar values, O is the least abundant element, followed by Ne. Relative to oxygen, 
iron has a significantly higher abundance, with Mg and Si a bit higher yet, and then Ni is 
much higher than the other elements. 
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3.3. Comparison of the RGS and EPIC Results 

We considered one-temperature models in both the central 1' (for both EPIC and RGS 
data) and more extended 4' diameter regions (for EPIC)/3.5' regions (for RGS), as well as 
two-temperature models (2T thereafter), shown in Tables 2-3. For the single temperature 
fits (Table 2), we give both the EPIC (MOS + PN) and the RGS results, and aside from 
the temperature and abundances, we list the degrees of freedom in the fit (dof), the reduced 
X 2 , and the luminosities from the power law component (Power) and from the thermal 
component in the 0.33 - 2.48 keV band (in units of 10 40 erg s" 1 ). The final column is the 
ratio of the thermal to the power-law luminosities. For the fit to the RGS data, we forced 
this ratio to be the same value as for the EPIC fit. In Table 3, we fit all the EPIC data 
and some good quality RGS data, but with a two-temperature thermal model, so there are 
two additional columns, for the second temperature, and for the luminosity from the second 
thermal component. The ratio of the luminosities is the sum of the two thermal components 
divided by the power-law component. There are two rows per instrument in each table, for 
the V (for both EPIC and RGS) and the 4' (for EPIC)/3.5' (for RGS) diameter regions. 

When fitting a model to the data, only photon statistics define the standard deviation 
per point. In addition, there are systematic uncertainties in the calibration of the instrument 
as well as instrumental artifacts, not all of which are known or excluded from spectral fits. 
There is no way of including accurately such effects in the spectral fitting, but we can identify 
the way in which the \ 2 will be affected. The resulting \ 2 will be larger than if there were no 
systematic effects. Furthermore, these systematic effects do not change as the source becomes 
brighter, but brighter sources have more photons and therefore lower uncertainties due to 
photon statistics. Therefore, systematic errors in valid spectral fits will lead to significantly 
higher reduced x 2 an d fits that formally would be rejected. This can lead to one rejecting 
spectral fits to all of the best data sets. Consequently, we have inspected every spectral fit 
and there are several that appear to be good fits but would formally be rejected. In Fig. [3] 
we show such an example of the best fit model (2T model) to the EPIC spectra of NGC 
4649, which has unacceptable reduced \ 2 values (1.350 for 1' region and 1.417 for 4' region) 
but shows no significant deviations between best-fit model and data. 

For the single-temperature fits within the central 1' (Table [2]), formally acceptable x 2 
values (not ruled out above the 99% confidence level) are found for three objects with fits 
to the EPIC data (NGC 720, NGC 4552, and NGC 5044) and two objects with fits to the 
RGS data (NGC 720 and IC 1459, although NGC 720 has the fewest RGS photons of any 
source). Of the sources with a large number of photons, which can be inferred from the 
number of degrees of freedom, the fit to NGC 4406, NGC 4649, and NGC 5044 appear to 
be good, while there are clear problems with the fits to NGC 1399, NGC 4636, and NGC 
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4472. The galaxy IC 1459 has a strong central point source that makes it difficult to isolate 
the thermal gas that we wish to study, so we exclude this from further discussion (but we 
gave the IT fitting results in Table 2 excluding central 13" region in radius). The fit using 
the 4' diameter region is worse than when using the central V region, so we emphasize those 
1' region results. 

In these one-temperature models, there is quite good agreement in the value of the 
temperature between the EPIC and RGS fits, and there is also good agreement in the Fe 
and O abundances between the EPIC and RGS fits within 90% error range except for NGC 
720, NGC 3923, NGC 4552, and IC 1459, which have very poor RGS spectra to constrain 
abundance parameters (Table [2]). Excluding the galaxies NGC 1399, NGC 4472, NGC 4636, 
and IC 1459 (see above), the range of the Fe and O abundances from the EPIC are 0.42-1.57 
and 0.26-0.55, respectively (relative to the Solar values). The median values for Fe and O are 
0.58 and 0.35 of the Solar values, while the median values from the RGS fits (also excluding 
NGC 720, NGC 3923, NGC 4552) are slightly higher, at 0.65 and 0.42. The ratio of the Fe 
to O abundance is around 2:1 in this data set. 

Of particular interest are three galaxies, NGC 4406, NGC 4649, and NGC 5044, where 
the RGS fits appear to be good and where the abundance uncertainties are not large (despite 
having formally unacceptable values of x 2 > discussed above). The galaxy NGC 4649 was 
discussed above and here we add that the RGS metallicities are in good agreement with 
those derived from the EPIC data. This is one of the few galaxies where the Fe metallicity 
is supersolar, although only modestly. The galaxy NGC 5044 has consistent abundances 
from the EPIC and RGS data, which show that the Fe abundance is about 0.74 of Solar 
and the O abundance is about 0.43 of Solar. The third galaxy, NGC 4406, has consistent O 
abundances (about 0.36 Solar) between the two XMM instruments, but the Fe abundance 
is higher in the EPIC observations. This is a complicated galaxy that is interacting strongly 



with its environment (IStickel et al.1 120031 ) . so perhaps it is not surprising that there would 



be disagreement for two instruments that sample slightly different spatial regions. 

When a second thermal component was introduced, the spectral fits generally improved 
and the metallicities usually increased. (We fit only six sources with good quality RGS 
spectra for a second thermal component added. The second temperature will be frozen to its 
EPIC value if it can not be constrained. We list all results but NGC 4649 in Table [3j which 
shows no convergence for the 2T model.) For the same eight galaxies used here (all except 
NGC 4636 and IC 1459) and for one-temperature fits, the Fe abundance rose by a median 
of 30% while the O abundance rose by 4%. For the EPIC fits to the data within the central 
V (Table ED, four of the galaxies had formally acceptable fits (NGC 720, NGC 3923, NGC 
4552, and NGC 5044), NGC 4406 is nearly acceptable, and the fits to the others appear to 
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be good, with the exception of NGC 4636. Of the good and formally acceptable fits (all 
except NGC 4636 and IC 1459), the range in the Fe abundance is 0.44-1.70 Solar and for O 
it is 0.27-0.76 Solar, with medians of 0.86 and 0.44 Solar. Even the two galaxies excluded 
from the analysis had similar values. The abundances for Mg and Si are indistinguishable 
from that of Fe, whereas in the one-temperature fits, they were systematically higher than 
the Fe abundance. In these fits, most of the power is from the gas phase, typically more 
than 90% in the 0.33 - 2.48 keV band. The RGS fits have Fe and O abundances of 0.33-1.50 
Solar and 0.31-1.22 Solar, with medians of 0.95 and 0.66 solar, which is consistant in Fe but 
higher in O. 

In Fig. H]- [9], we show zoomed-in RGS spectra with best fit model for the six most bright 
galaxies with marked individual emission lines. Note that, overall the fits are fairly good, 
but there are some residuals at Fe XVII and Fe XVIII at 15 A, 16 A, or 17 A. Those residuals 
are not due to improper convolution with large image size, since usually the isolated O VIII 
line at 19 A is fitted fairly well. They might be due to calibration issues or non-accurate 
APEC model at those lines. 



One important object for comparison is NGC 4636, where Ku et al.l (120021 ) fit a model 
to the RGS data within 2' that included a temperature gradient and optical depth effects 
for Fe XVII, which necessitated the use of a ft model. Their ab undances for Fe and O ar e 



1.29 ± 0.06 and 0.63 ± 0.06, relative to the Solar abundances of iGrevesse fe Sauvall (119981 ). 
Their ratio of Fe to O is 2 in these units. From our fitting of the RGS data with a single 
temperature and no opacity correction, we obtain Fe and O abundances that are about 30% 
lower, but with the same abundance ratio. Our so mewhat lower a bundance for Fe might 



be due to the absence of optical depth corrections. Ku et al.l (120021 ) found that the optical 
depth of Fe XVII line at 15.0 A is greater than unity for the densities and temperatures in 
the core of NGC4636 ISM, while the Fe XVII blend at 17.1 A is negligible. Many scatterings 
of the 15.0 A photons flatten its profile comparing to the others, which would cause large 
uncertainty on Fe abundance. To correct the optical depth effect at 15.0 A, we fit again 
the 3.5' apeture RGS spectra of those very bright galaxies (NGC1399, NGC 4406, NGC 
4472, NGC4636, NGC 4649, and NGC 5044), excluding the problematic 15.0 A line (i.e. 
excluding 14.75 - 15.4 A ) as well as freeing the redshift parameter to account any residual 
line shift. It turns out that only NGC 4636 shows large changes in both parameter values 
(18% increase in Fe, and 33% increase in O) and in reduced xl (decrease by 23% ). Other 
galaxies only show at most 4% variances in Fe, O, and reduced Xr- So optical depth effect at 
15.0 A line is negligible for all the galaxies except for NGC 4636. After corrected this effect 
for NGC 4636, we still have slightly lower abundances of Fe (1.07±g;^) and O (0.57lg;{$) 
from IT fit, but have consistent abundances of Fe (1.421q' 17 ) an d O (0.58 + ^) from the 



2T fit, comparing to Fe (1.29 ± 0.06) and O (0.63 ± 0.06) from Ku et all (12002( 1 . These 
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abundances are slightly higher than what XMM EPIC data and ASCA data reveal. iBuote 
Jl999h obtained abundances of Fe and with ASCA of 1.08±°;g and 0.49±Jj;?jj, respectively, 
which are consistant with our 2T fit for XMM EPIC data at 4' diameter region, 1.06 ± 0.03 
for Fe and 0.51 ± 0.03 for O. The reason that these abundances from CCD spectral analysis 
are systematically lower (25% lower for Fe and 12% lower for O) than the XMM RGS results 
is not due to the uncorrected optical depth effect for Fe XVII line at 15 A, since we refitted 
the XMM EPIC 4' diameter region excluding this line (i.e. excluding 0.79 - 0.88 kev), and 
found only 4% increase in Fe and 2% increase in O. The reason might be that the RGS 
and EPIC cameras cover different regions o f this galaxy, which shows complex struct ure and 
events that have been described by others (jjones et al. 2002 ; O 'Sullivan et al. 2005 ). 



3.4. Comparison with other studies of the metal abundances in early-type 

galaxies 

We show below the detailed comparison with other studies of the metal abundances 
in early-type galaxies for individual galaxy. Note that our fitted temp eratures agree well 
with t he published results, and abundances from literatures are scaled to iGrevesse fe Sauval 
jl998h . 



NGC 720 - Our 4' diameter aperture EPIC results for NGC 720 (Fe(l.lO ^), O(0.45j£ff 



Mg(0 : 76 



h0.73\ 
-0.32/ ' 



and Si(0.97_ ; 45 )) are consistent with Chandra results reported by lHumphrey fe Buote 



(120061 ) (Fe(0.71 



!-0.40\ 
-0.2lJ) 



O(0.16 ± 0.16), and Mg ( 0.90 ± 0.37)), and are also consistent with the 



180 ks Suzaku observation by ITawara et al.l (120081 ) within 90% error range, who obtained 



abundances of Fe(0.73ig;££), O(0.47ig;}£), Mg(0.50lg;&), and Si(0.54l°; 



NGC 1399 - T his is an imp ortant galaxy with supersolar abundances reported by pre- 
vious observations. iBuotd ( 19991 ) found an Fe abundance of 2.40l 9 o within a 10' diameter 
aperture by ASCA. Later, he reported an XMM result fo r this galaxy that the Fe abundance 
in the central region (about 10' in diameter) is 1.54-2.04 ()B_uoteJl2002J). Our net 120 ks XMM 
observation also shows a consistent but more constrained Fe abundance, 1.32 ± 0.05, within 



4' diameter region. This Fe abundance is higher than what iHumphrev fe Buotd (12006!) ob- 
tained with Chandra data, 1.06 ±0.09. This is because IHumphrev fe Buotd (120061 ) obtained 
an emission-weighted average Fe abundance up to 20' in diameter, which tends to be lower 
becau se the Fe abundance decr eases to subsolar at large radius (20' in diameter) of NGC 
1399 IHumphrev fe Buotd lioo3 l. 



NGC 3923 - IHumphrev fe Buotd (120061 ) reported a poorly constrained abundance of 
Fe for 1.03(> 0.24) due to very few counts of photons. We obtained a lower but more 
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constrained Fe abundance (0.72loig)- Since we have many more counts of photons (total 
counts of 3.3 * 10 4 ) for this galaxy and obtained a good fit (xf — 1-024), we are confident 
with our result. 



NGC 4406 - ASCA observation shows an Fe abundance of 0.77±g;|f flBuote et alJll998h . 



while Chandra observation shows Fe and O abundances of 1.10i ; 4 Y and (UOioit, 

respectively 



(lAtheyi 120071 ). Our results are consistent with these two within error range, which are Fe of 



0.87±°;°;j and O of 0.51±g;g|. 



NGC 4472 - iBuotd 
for Fe abundance, while 



fl!999h reported an ASCA result for this galaxy to be 2.961^ 



Humphrey fe Buotd (120061 ) reported an emission-weighted average 



abundances for Fe and O of 1.25l ; 3 6 an d 0.48lo!i8; respectively. Our XMM result gives 
consistent but more constrained Fe and O abundances of 1.75^ao6 anc ^ O-^-aoL respectively 



NGC4552 - [Humphrey fe Buotd (12006!) measured abundances of Fe and O from Chandra 



data of 0.63j£o£ and 0.08 ± 0.06, while lAtheyi (120071 ) reported abundances also measured 



from Chandra data up to 2.5' diameter region to be 0.73±g;g? for Fe and 0.29lg;?£ for O. We 

0.05 > respectively, 



-0.67 ■ L C allyj - w - ic, -0.18 

obtained better constrained abundances of Fe and O of 0.49±0.08 and 0.25 +0 



which are consistent with lAtheyi (120071 ) 's Fe and O abundances w ithin 90% error range, 
but is 68% higher in O abundance than Humphrey fe Buotd (j2006l ) , s result. We attribute 
this discrepency to be not having enough counts in Chandra's spectra to constrain the O 
abundance. 

Extra comment for NGC4552: This system is known to be suffering fairly strong ram- 
pressure stripping and the extent of its gas is quite well-defined from the Chandra images 
(IMachacek et al.l 120061 ) . The 1' diameter region is quite a good match to the gas, but the 
4' aperture might contain a large fraction of cluster and mixed cluster /galaxy gas. In our 
two-phase model fits in Table 3, however, the abundances are quite consistent between 1' 
and 4' regions. Therefore, the cluster emission will not significantly affect the measurement 
of the galaxy hot gas abundances in NGC4552 even for 4' apeture. The reason may be that 
even for the 4' region, the local hot gas in NGC4552 is still the dominant thermal component 
comparing to other thermal component (like the cluster emission), which contributes 63% of 
the total thermal luminosity. 

NGC 4636 - See detailed discussion in Section 13.31 

NGC 4649 - See detailed discussion in Section 13.21 



NGC 5044 - iTamura et al.l (120031 ) obtained abundances of Fe and O within the central 
10-20 kpc region (2' width) of NGC 5044 from XMM-Newton RGS data. Their abundances 
for Fe and O are 0.81 ± 0.07 and 0.31 ± 0.13, respectively, which is consistent in O with but 
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higher in Fe than our RGS 1-T model fitting results for V region, Fe(0.65l[j;{$), O(0.42^°), 
but ar e consi stent with our RGS 2-T model fitting results of Fe(0.77+°;^) and 0(0.48+°;}?). 



Atheyl (120071 ) used Chandra for annualar analysis for this galaxy, and found abandances of 
Fe, O, Mg, and Si within 4' diameter in the ranges of 0.92-1.53, 0.18-0.96, 1.01-1.77, and 
0.95-1.93, respectively. These abundances gradients are also conformed by our study with 
two aperture analysis of EPIC data. We find abundances gradients of Fe, O, Mg, and Si 
betweem smaller and larger apertures in the ranges of 0.86-1.15, 0.41-0.49, 0.84-1.14, and 
0.86-1.01, respectively. Our a yerage Fe abundan ce (1.15 ± 0.05) within 4' diameter is also 
consistent with the result by iBuote et al.l (120031 ) . who reported the XMM observation of 
emission-weighted average abundance of Fe within 10' diameter to be 1.09 ± 0.04. 

IC 1459 - This is the only galaxy in o ur sample that shows different spectral charactis- 
tics due to a strong central point source (jFabbiano et al.ll2003l ). Its power-law component 
is extremely large compared to the thermal component, making it hard to constrain the 
abundances in the thermal emission. We fit its spectra excluding the central 13" region 
in radius, which is mainly a power-law c omponent with a pho ton index V of 1.88 and an 



excessive absorption rin of 2.9 x 10 cm (jFabbiano et al.ll2003l ). We then fit the remainder 



(mainly thermal emission) with a 1-T thermal component plus a power-law (phont index 
T fixed to 1.56) with only Fe allowed to vary (other elements are tied to Fe). This gives 
a better constrain on the Fe abundance of 0.39+Qog, which is consistent with the Chandra 



result of 0.74to 6 u (jAthevi 120070 and the ASCA result of 0.21+°,;^ feuote et al.lll998fl . within 
their error ranges. Longer observation is needed to constrain other element abundances. To 
test if the central point source contributes any photons outside the excluded 13 " region, we 
refitted IC1459 EPIC 1' region excluding the central 13" core with photon index free to fit. 
The fitted photon index is 1. 60^05, only increased by 2.6% from 1.56, and the temperature 
and Fe abundance change not much, by 0.5% and 3.7%, respectively. For 4' region, the 
results are essentially not changed. So we conclude that the central point source in IC1459 
doesn't contribute many photons outside the excluded 13" region. 



3.5. Discussion on the Assumptions of the Spectral Fitting Method 

In our spectral fitting method, we made two simplifying assumptions in the choice of 
spectral models - freezing the slope of the power law component and linking the abundances 
of the two-temperature models. We did the following tests to justify these assumptions. 

Effect on the slope of the power law component: To check the effect on pho- 
ton index, we test two cases (with 2T model plus a power law), NGC5044 and NGC3923, 
according to their different fractions of the total luminosity for the power-law component. 



- 15 - 



During the fitting, the photon index is free to fit. 

For NGC5044, the power-law component only contributes 3.4% to the total luminosity 
in the energy band of 0.33 kev to 2.48 kev in its central 1' diameter region. The fitted 
parameters don't change much. The result shows a best-fit photon index of l.SO+g^g, which 
is consistent with 1.56 within its 90% error range. The two temperatures vary by only 0.3%, 
while the abundances vary at most 5.8%, which are within 90% confidence error range of 
previous result when fixed photon index to 1.56. 

For NGC3923, the power-law component contributes 24.6% to the total luminosity in 
the energy band of 0.33 kev to 2.48 kev in its central 4' diameter region. The fitted photon 
index is 1.73+0^, which is consistent with 1.56 within its 90% error range. There are no 
significant changes in reduced chi-sqaure value as well as the two temperature values, only 
by 0.4%. Abundance values, however, show large increase by at least 50%. Fe changes most, 
increased by 70.8%, from 0.72+°;^ to 1.23+ ^44, but is still consistent within their error range. 

So we conclude that the value of photon index doesn't significantly affect the abundance 
values for hot gas dominent galaxies, and freezing photon index to 1.56 can better constrain 
the abundance values. But it is the main error source for galaxies where power-law com- 
ponent is comparable to the hot gas component. In our sample, most sources only have a 
few percent of the total luminosity for the power-law component, and thus are not greatly 
affected by the value of photon index. Only four galaxies (NGC720, NGC3923, NGC4552, 
and IC1459) have power-law component with luminosity greater than 15% of the total lumi- 
nosity, thus their abundances are more greatly affected by the value of photon index. Longer 
observations for these sources are needed to better constrain the abundances. 

Effect on linking the abundances of the 2T models: There are five galaxies 
(NGC1399, NGC4406, NGC4472, NGC4649, and NGC5044) in our sample which show the 
possible signs of overlying cluster emission. This can be seen from the 2T model fitting 
temperatures (see Table 3), where the galaxy hot gas temperatures are around 0.63-0.81 
kev, while the cluster emission temperatures are around 0.94-1.72 kev. 

We chose NGC1399 4' region to test the effect of linking two abundance sets, because 
it has the most photon counts and a comparable cluster emission component to the local 
galaxy hot gas (about 46% of the total thermal luminosity). During the fit, we unlinked the 
two abundance sets. Each abundance set then has free parameters of Fe, O, Mg, Si, S, Ne, 
and Ni. Other elements are tied to Fe. The fitting results show no change for the reduced 
chi-sqare value, but the 90% confidence range for each abundance parameter becomes larger. 
The best-fit parameters for the lower temperature component are T:0.81llo;oo4, F e: l-89^a57> 
0:1.29^0:69, while values for the higher temperature componet are T:1M5±%^, Fe:1 - 20 -ai3, 
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and 0:0.101 

o!io- These two sets of abundances do not show much difference within their 
90% error bars. They are also consistent with the best-fit parameters when linking the 
abundances of the 2T models. 

So in our sample, even for the region that has the most contamination from the cluster 
emission, linking the two abundance sets would not significantly affect the hot gas abun- 
dances, and it can better constrain the abundances by using fewer free parameters. Direct 
measurement of abundances for different thermal components are not good for low quality 
spectra, and deprojected method should be used. 



3.6. Discussion on the Effect of RGS Background 

For the 1' RGS region, the local galaxy hot gas emission is always the dominant compo- 
nent. Thus it makes little difference, whether using an RGS background template or using 
one that contains the underlying cluster emission but free of lcoal galaxy hot gas emission. 

For the 3.5' RGS region, however, the background could have some effect, particularly for 
those with the brightest surrounding group/cluster emission. In our sample, fainter systems 
like NGC4552, NGC720, NGC3923, and IC1459 have too poor RGS spectra to constrain 
their abundances, even for the 3.5' region. For brighter systems, however, it is very hard to 
get a local background that contains only the underlying cluster emission but free of galaxy 
hot gas emission. This is because their extended hot gas emission regions are usually larger 
than 5' in diameter (checked with XMM MOS images), but the maximum field of view of 
RGS in the cross-dispersion direction is only 5'. The only way to obtain this local background 
is using seperate observation pointing off the target and on the surrounding cluster emission, 
which is not available for all the galaxies in our sample. Thus using local contaminated 
background will overestimate the background, while using the RGS background template will 
underestimate the underlying cluster emission. This is the main error source when determine 
the abundances for the RGS 3.5' region, especially when the cluster hot gas component is 
comparitable to the galaxy hot gas component, which can be seen by the luminosity ratio of 
the two thermal components obtained from the EPIC fits. 



4. Discussion and Conclusions 



There have been several previous determinatio ns of metal abundances in early-type 
galaxies, som e of which have obtained similar results (jBuotelll999l ; iHumphrey &: Buotdl2006l ; 
Atheyl 120071 ). One important difference in this work is that, for a sample of 10 galaxies, we 
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include the RGS data, in which one can clearly see the strong individual lines from Fe, O, and 
other species, whereas they are thoroughly blended in the XMM EPIC and in the Chandra 
ACIS data. The identification of these ions not only gives us greater confidence in the 
abundance determinations, but they often provide very strict constraints on the maximum 
allowable abundances. 

Of the models considered, the one yielding the lowest x 2 has a power-law to represent 
unresolved X-ray binaries, two thermal components with one set of variable abundances for 
several individual elements (Fe, O, Mg, and Si), and Galactic absorption. This model led 
to formally acceptable or good fits for most of the galaxies except for NGC 4636 and IC 
1459. In IC 1459, there is a strong central point source that compromised spectral fits to the 
fainter emission from the hot gas. In NGC 4636, a two-temperature model is inadequate, 



possibly due to the structu re and events that have been described by others (IJones et al. 



2002l ; IO , Sullivan et al Jl2005h . Excluding NGC 4636 and IC 1459, the remaining eight galaxies 
have median abundances (2 5% and 75% quartile poin ts are given in parenthesis for Fe and 
O) relative to Solar values JGrevesse fe Sauvall Il99sf ) of 0.86 (0.64-1.70), 0.44 (0.29-0.60), 
0.81, and 0.79 for Fe, O, Mg, and Si. The abundances for Mg and Si are not distinuishable 
from that of Fe. Three of the galaxies have mildly supersolar Fe abundances, at about 1.5 
of the Solar value, and we note that the three all lie at the center of their group or cluster 
(NGC 1399 in the Fornax cluster; NGC 4472 is the dominant galaxy in the southern part 
of the Virgo cluster; NGC 4649 is the brightest member of a group on the outskirts of the 
Virgo cluster). While this may be a trend, we note that the galaxy NGC 5044 is also the 
brightest member of its group yet it has abundances that are close to the median values. 
The galaxy that has the lowest Fe abundance is NGC 4552, the optically least luminous 
galaxy of this group (and relatively lower in X-ray luminosity). As we will discuss in a 
future paper, relatively low abundances are found frequ ently in galaxies with lower X-ray 
and optical luminosities, as lO'Sullivan fe Ponmanl (120041 ) has found. 



One of the issues that we examined was the difference between abundances when a sec- 



ond temperature c omponent is added, sometimes referred to as the Fe-bias (jTrinchieri et al. 
19941 ; iBuotd |2000| ). The need for a second temperature component is probably due to a 



temperature gradient within the galaxy. Upon adding the second thermal component, the 
increase in Fe was about 30% (quartile values of 7-38%) and for O it was 4% (quartile values 
of -3%~16%). There are two galaxies in which there was essentially no change in the abun- 
dances, NGC 4649 and N GC 5044. F r om radial temperatu re profiles, NGC 4649 is known 
to be nearly isothermal (jAtheyl 120071 ; iRandall et al.ll2006l ). which explains why adding a 
second temperature component was of little consequence. For NGC 5044, there was always 
a dominant thermal component and the luminosity from an additional thermal component 
was relatively small in the best-fit solutions. 
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Because NGC 4649 can be adequately fit with a single-temperature component and the 
temperatures do not change with an additional thermal component, we used this source to 
compare results between spectral fits for the three XMM instruments (EPIC MOS, EPIC 
PN, and the RGS) and for the Chandra ACIS-S. The derived temperatures were the same to 
about 1%, and there was excellent agreement between the three XMM-Newton instruments 
for the seven fitted abundances. The abundances derived from the Chandra ACIS-S were 
also in good agreement for Fe, Mg, Si, and Ni, but gave lower values than the XMM-Newton 
data for O, Ne, and S. Overall, Fe, Mg, and Si had similar abundances that were 1.5 times 
the Solar value, S is consistent with that value, and Ni may be supersolar, although the 
constraints are poor. Oxygen, the lightest metal in the group, has an abundance of about 
0.6 Solar and Ne, the next lightest metal, has an abundance between O and Fe, at about 0.7 
Solar. 

Metals in hot gas of early type galaxies are thought to come from stellar mass loss, 
along with supernovae. Since supernovae mainly contribute heavy elements to the hot gas, 
we should expect supersolar metalicity in hot gas. In Fig. [TU1 we compared Fe abundances 
found in hot gas and those from optical stud ies (only six galaxies in ou r sa mple have op- 



tical F e abundances reported in literature by Humphrey fe Buotd (120061 ) and iTrager et al. 



( 12000( 1 ). Only three galaxies (NGC 1399, NGC 4472, and NGC 4649) show modest super- 
solar Fe abundances in hot gas, while the other three (NGC 720, NGC 3923, and NGC 
4552) are slightly subsolar. Supernovae were expected to enhance the gas phase abundances 
considerably, which is not observed in our sample. 

Element abundance ratios can constrain the relative contributions of SN la or SN II 
to the metallicity. In our sample, the median abundance ratios are 0.51, 0.94, and 0.92 for 
O/Fe, Mg/Fe, and Si/Fe, respectively. These give us the fraction of contribution from Type 
la supernovae of 6 4%-85% when u sing l inear combination of SN la yeilds from W7 model and 
SN II yields from iNomoto et al.l (119971 ). which indicates that most of the metal enrichment 
is from Type la supernovae. 

One important result is t he low O/Mg rati o, which we derived a valu e of 0.54, which 
is also confirmed to be low by lAtheyl (120071 ) and iHumphrey fe Buotd ( 120061 ). Theoretically, 
One should expect similar abundance of O to Mg, since SN la produces similar amount of O 
and Mg while SN II produce slightly more O than Mg. The chemical studies of Galactic bulge 
stars, which are also old star population similar to stars in elliptical galaxies but are mainly 
enriched by Type II supernovae, also show a si milar low O/Mg ratio at high Fe abundance 
( jFulbright et al.l 12004 ; iMinniti fe Zoccalil 120071 ) . This low O abundance in the hot gas as 
well as in bulge stars might be the evidence of an overestimate of O yield by SN II models, 
which do not consider significant mass loss at the late stage of massive progenitor stars. 
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Since strong stellar wind from massive progenitor stars might reduce or even completely 
blow away their He-rich envelopes while still preserve the C-0 layers, thus after Type II 
supernove explosions, only a elements (like Mg) by Carbon burning are synthesized whil e 
elements like O through hydrostatic He burning are greatly reduced ( IFulbright et all 12004 ) . 
The obser ved Mg-rich SN re mnent of N49B without accompanying Ne or O can support this 
argument (IPark et al.1 120031 ). 



We would like to thank Renato Dupke for his comments and suggestions. We gratefully 
acknowledge financial support for this research, which was provided by NASA. 
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Fig. 1. — XMM-Newton full-range RGS spectrum of NGC4472 with marked emission lines 
from 7 A to 26 A band. 
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Fig. 2. — The temperature and abundances from the XMM-Newton instruments and the 
Chandra ACIS-S for NGC 4649. A 1' diameter region was used, along with a model that 
includes Galactic absorption, power-law emission from binary sources, and a one-temperature 
thermal model where abundances are fitted for the elements shown. The results from the 
XMM-Newton instruments are in agreement with each other and with the Chandra results, 
with the exception of O and possibly S. Elements heavier than Ne are modestly supersolar, 
while the two lighter metals, O and Ne, have abundances about a factor of two lower. 
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Fig. 3. — Example of NGC4649 for XMM-Newton EPIC spectral fitting with marked emis- 
sion lines. The upper panel is for 1' aperture region while the lower panel is for 4' aperture 
region. (Green) :pn camera; (red): MOS1 camera; (black): MOS2 camera. Note: the emis- 
sion lines from Fe L complex and O VIII are completely unresolved while emission lines from 
Mg, Si, and S can be indentified fairly well in CCD spectra. 
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Fig. 4.— The 10-20 A region of the XMM-Newton RGS1 and RGS2 spectra for NGC 1399 
with the best-fit models. The upper pannel is for V aperture while the bottom one is for 3.5' 
aperture, with marked resolved emission lines mostly from the Fe L complex and O VIII. 
Note that the wide gap between 10.5 A and 13.8 A is due to failure of RGS1 CCD7, while 
the narrower gaps at 11.4 A, 13.1 A, 14.6 A, 15.3 A, 16.5 A, 17.1 A, and 18.9 A are due to 
gaps between RGS CCDs. 
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Fig. 10- — Comparison of Fe abundances from hot gas (our results) and from optical studies 
of stars ( Humphrey fe Buote 2006] ; Trager et al. 2000 ). Symbol X: mean stellar Fe abun- 



dances from 



Trager et al. 



Humphrey fe Buotd (120061 ); symbol diamond: mean stellar Fe abundance from 
(12000 ). The solid line represents where the hot gas and the stars would have 



the same Fe abundance (it is not a fit to the data). 



Table 1. Properties and Observational Information of Surveyed Galaxies 



Galax T c a S ?" b D ° Nh A r "° Ob Id' Nct ox P osurc timc ( ks ) E 

a axy ype (mag) (Mpc) (10 20 cm -2 ) (arcsoc) '"' MOS1 MOS2 PN RGS1 RGS2 



NGC720 


E5 


11.13 


27 


67 


1 


.55 


39.87 


0112300101 


29.6 


30.0 


19 


9 


46.3 


45.3 




NGC1399 


cD;Elpcc 


10.44 


19. 


.95 


1 


31 


42.55 


0400620101 


121.5 


120.5 


73. 


9 


118.2 


118.1 




NGC1399 
















319 h 














55.9 


NGC3923 


E4-5 


10.62 


22 


.91 


6 


.29 


53.35 


0027340101 


38.8 


38.7 


29 


8 


43.9 


42.6 




NGC4406 


S0(3)/E3 


9.74 


17. 


14 


2 


.58 


89.64 


0108260201 


77.9 


78.9 


47 


8 


83.5 


81.1 




NGC4472 


E2/S0 


9.33 


16 


.29 


1 


.65 


104.40 


0200130101 


82.5 


82.7 


72 


8 


101.5 


101.5 




NGC4472 
















321 h 














32.5 


NGC4552 


E 


10.57 


15 


.35 


2 


.56 


48.89 


0141570101 


27.8 


31.2 


18 


5 


42.9 


42.8 




NGC4636 


E/S0.1 


10.43 


14 


66 


1 


.83 


100.08 


0111190701 


59.2 


59.3 


51 


1 


62.9 


61.4 




NGC4649 


E2 


9.70 


16 


.83 


2 


.13 


73.73 


0021540201 


50.6 


50.6 


42 


2 


53.1 


51.6 




NGC4649 
















785 ft 














22.9 


NGC5044 


EO 


11.67 


31 


.19 


5 


.03 


82.23 


0037950101 


22.6 


22.7 


17 





23.6 


22.8 




IC1459 


E3 


10.83 


2!) 


24 


1 


.19 


38.61 


0135980201 


29.3 


29.3 


25 


2 


31.7 


30.8 





a Thc galaxy type was take from NED. 

b Total B band magnitude from RC3 (Do Vaucoulcurs ct al. (1991)). 

c Distances in Mpc, measured by surface brightness fluctuation method ( Tonry ct al. (20(31)). 
d The Galactic H I column density, taken from the dust map by Dickey & Lockman (1990). 
c Effcctivc, blue-half light radius in arcsecond derived from R.C3 |Dc Vaucoulcurs ct al. (1991)). 

XMM-Newton observation ID. h -.Chandra ObsID. 
g Nct exposure timc after filtering background flares. 
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Table 2. One Phase Model Fitting for Chandra ACIS-S, XMM EPIC, and RGS spectra. 



Galaxy 



Camcra/Sizc a 



T(keV) b 



Mg c 



dof d x^ c Power f Vapec s Vapec/Power h 



NGC720 



EPIC 1' 
EPIC 4' 
RGS 1' 
RGS 3.5' 



0.5607 



,,+ 0.015 
-0.008 
QQfi + 0.013 
U.dMb_ Q 021 

0.369l°;J^ 
0.327 + " " 



0.51 
1.52 



+ 0.06 
0.05 
+ 1.27 
0.41 
1+3.46 



0.30Z 
0.49 4 



4.10^ 



,+0.17 
-0.13 
, + 0.54 
0.16 

010 ±o.ou 
0.31+2?? 



0.40 
1.35 



+0.19 
0.16 
+ 1.54 
0.46 



0.35 
2.48 



+0.30 
0.25 

+ 2.52 
1.09 



0.54 
2.21 



+ 0.30 
0.29 
+ 0.15 
0.15 



5+1.02 
J -0.98 



198 1.102 

500 1.274 

62 1.313 

143 1.303 



0.54 
1.61 

0.92 
2.06 



u.ouy_ 003 
n „, , + 0.002 

O qc n + n - 002 

O R1Q + 0.011 

0.813_ „ 14 
R + 0.008 
r ' -0.009 



i ,., + 0.10 

i - w -0.09 
, ,- + 0.04 
1 -*°-0.0i 
„ fi - + 0.02 
uoo _0.02 
, ,, + 0.62 
lc " i -0.37 
+ 0.10 
0.06 



77+11-14 
u ' ' ' -0.12 

54+ 07 
u -°^_0.07 

50+ 04 
u ' ou -0.03 

1 20+ 029 
lzu -0.36 
n ss+0-10 

0.07 



1 49+U-^ 
^ ^ -0.19 

1 65 + 011 

89+° 04 
u -° y _0.04 



„ + 0.20 
-0.18 
, + 0.09 
"-0.09 
, + 0.04 
u -0.04 



37+U-36 
u<s '-0.34 

1 97+0-13 
L - y < -0.12 
„ ,,+0.08 
u -^-0.04 



°- 76 -0.28 

1 13 + 011 
1 - lcl -0.12 

89+ 03 



,+ 1.02 

"-0.87 



7+0.11 
' -0.11 



1.73 
4.19 
2.81 
5.37 
7.98 
7.68 
19.94 
14.12 
26.88 
3.48 
5.20 
3.76 
4.66 
0.08 
0.44 
0.29 
0.80 
4.15 
4.36 
10.50 
8.37 
13.56 
0.18 
0.22 
0.18 
0.25 
4.83 
14.89 
6.93 
14.57 
7.19 
14.85 
5.68 
5.66 
9.63 
6.93 
9.41 
49.05 
249.84 
123.97 
293.81 
0.31 
1.36 
1.95 
2.82 



3.22 
2.61 
3.05 
2.61 
11.28 
12.18 
5.77 
16.46 
5.91 
5.68 
2.80 
5.86 
3.24 
8.26 
12.07 
8.32 
13.29 
14.11 
11.01 
7.59 
15.01 
9.06 
2.64 
1.69 
2.59 

I. 67 
34.70 
27.92 
34.92 
27.91 
36.22 
28.48 

II. 91 
9.73 
5.50 
9.73 
5.67 
26.87 
28.39 
26.95 
28.65 
0.35 
0.69 
0.97 
0.99 



ACIS-S l' 
EPIC 1' 
EPIC 4' 
RGS 1' 

RGS 3.5' 



0.996 ■ 



0.68 



47 +0.016 

n ,07+0.016 
0.387_„ Q14 

446 + " 050 

432+ 057 
-0.067 
64 7+lUHIY 
U - D *' -0.007 
, ,+ 0.003 
°- 716 -0.003 

°- 626 -0.018 
623+ 011 



4 9 +l'-l )!J 
u '* a -0.04 

63+ 017 
U ' OJ -0.12 

5 00+ 00 

5 00+ 00 
J, " u -3.33 



°' 26 -0.07 
21+ 06 

2.06 



044+^1 0.98±g|g 
0.45±!H 4 0.63+1?* 0.27±!H 4 



143 3.520 
775 2.393 
1367 2.939 
1011 1.190 
1451 1.323 



0.71 
0.63 
3.46 
0.86 
4.55 



' , ' 11 -0.78 
, 7 ,+ 1.07 
£ - ' J -1.08 



NGC3923 



EPIC 1' 
EPIC 4' 
RGS 1' 
RGS 3.5' 



+ 1.33 
1.06 



-0.10 
, + 0.22 
"-0.16 



-0.19 
, + 0.48 
u -0.35 



0.14 
+0.10 
-0.10 



3.01 



+ 1.! 



315 1.346 

588 1.168 

91 1.421 

184 1.476 



0.61 
1.86 
0.64 
1.44 



+ U.04 

0.04 
+ 0.03 
0.03 
+ 0.06 
0.04 
+ 0.03 
0.02 

— -rn-nr 

J-'OO-O.IS 
, ,,+0.06 
^•^-0.03 
+ 0.03 
0.03 
+ 0.17 
0.11 
+ 0.29 
0.15 

u -*^-0.05 



65+ U 11 
u - DO -o.io 

65+° 05 
U - DO -0.05 

15+ - 20 
u - 1;> _0.05 

0.16 ' 



+ 0.13 
0.12 

+ 0.05 
0.05 



+ 0.19 

0.20 
+0.10 
0.10 



+ 0.46 

0.52 
+ 0.24 
0.22 



EPIC 1' 
EPIC 4' 
RGS l' 
RGS 3.5' 



0.64 
0.65 
0.34 
0.29 



40 T 
u ' w -0.09 

49+ 05 
u ' w -0.05 

32+ 11 



0.69 
0.62 



0.83 
0.82 



, + 0.09 
'-0.10 



1.79 
2.33 



0.25 



724+ lu,1!) 
u. '^ 4 _o.008 

O 748 + " 003 
U - 748 -0.003 
„ ,-,, + 0.002 
O.S13_ 002 

„ 7 ,„+o.oio 

O o, O+0-003 

u.»i/_ n , )06 



0.08 
+ 0.05 
0.04 
+ 0.29 
0.23 
+ 0.07 
-0.06 
+ 0.04 
0.04 
+ 0.14 
0.05 
+ 0.24 
0.14 
+ 0.09 
0.08 
+ 0.08 
0.06 
+ 0.39 
4.44 
+ 1.59 
-0.96 



+ 0.11 
0-06 



462 1.273 

919 1.409 

304 1.294 

774 1.230 



0.01 
0.04 
0.03 
0.06 



, + 0.45 

-0.34 
x 17 +0.09 



, O7+0-41 
i - S '-0.32 
1 19 + 009 
- L - ly -0.08 
, ,, + 0.06 

-•-■•"-o.os 



r + l-H4 

J -0.86 
+0.15 
0.14 
+ 0.09 
0.08 



+ 0.58 

0.49 
+ 0.13 
0.13 
+ 0.08 
0.09 



, ,r- + i K 3 
• 3 ' 50 _1.59 

3 93+ 32 



ACIS-S l' 
EPIC l' 
EPIC 4' 
RGS l' 

RGS 3.5' 



1.05: 
0.71 4 
1.04 4 



0.71 
0.54 
0.51 
0.59 
0.91 



1.99"* 



1.55 
0.53 
0.81 



«.608i!!;ii;ii 

o-386t°;°i? 
o.557t° ;°of 4 

0-3 7 6lg:o 4 4 9 o 



, 08 
+ 0.07 
0.06 
+ 0.31 
0.18 
+ 0.48 
0-27 



0.83 
2.81 



1.00 
0.99 
1.18 



123 1.595 

720 1.771 

1162 5.235 

536 1.478 

1015 1.231 



0.29 
0.40 
1.38 
0.56 
1.50 



46+" 11 44+" 13 

U - 4D _0.10 u - 4 -0.12 

; +0.33 , 4c ,+0.72 

0.23 - L '*°-0.53 



NGC4552 



EPIC l' 
EPIC 4' 
RGS l' 
RGS 3.5' 



0.66 
4.32 
5.00 



„ -..,- + 0.006 

631+ 002 
u - oal _0 .001 

622+" 1,09 
639+ 006 
n eo Q +0 0<>9 

°- 589 -(),008 
°- 62j -0.006 



+ 0.20 
0.12 
+ 0.68 
3.33 
+ 0.00 
-3.24 



0.29 
0.22 
4.61 
2.59 



o.ycr 



, + 0.07 
°-0.07 



414 1.063 

628 1.295 

89 1.640 

190 1.469 



0.07 
0.13 
0.07 
0.15 



+ 0.01 

0.02 
+ 0.03 
0.03 
+ 0.10 
-0.10 
+ 0.10 
0.07 
+ 0.12 
0.0. f 



+ 0.04 

0.03 
+ 0.03 
0.03 
+ 0.08 
0.07 

43 + 07 
u ' M -0.06 
, + 0.04 
0.08 



+ 0.04 

0.04 
+ 0.03 
0.03 
+ 0.18 
0.16 
+ 0.16 
0.13 
+ 0.14 
0.21 
+ 0.12 
0.17 
+ 0.29 
0.28 
+ 0.18 
0.16 
+ 0.12 
0.11 



+ 0.06 

0.06 
+ 0.06 
0.06 



+ 0.31 

0.25 
+ 0.19 
0.20 



EPIC 1' 
EPIC 4' 
RGS l' 
RGS 3.5' 
RGS corrected 1 1 
RGS corrected 3.5' 



0.62 
0.92 
0.83 
0.91 
1.00 



0.41] 
0.58j 
0.38 4 



o, r ,:v 



763+ UU14 
77R+0-003 
u. / "J_o,003 

O 784 + 0- 002 
U - 784 -0.002 
n 7CQ + 0-010 
U. 'o«_ oio 

776+0 014 
u ' ' '-0.009 



1 07+ 011 
lu '-0.02 



+ 0.21 

0.19 
+ 0.12 
0.11 
+ 0.08 
0.08 



„ - 7 +0.03 
u - J '-0.07 



0.40 
0.69 
0.45 
0.48 
0.65 
0.63 



2.35 
1.06 



506 
838 
438 
893 
389 
810 



2.061 
2.167 
2.020 
1.700 
1.342 
1.307 



0.14 
0.53 
0.20 
0.52 
0.20 
0.52 



ulo -0.13 

55+ 11 
u ' oo -0.10 
„ c-, + 0.08 
°' 5S -0.07 
0.67+ ' 33 



! 47+0-iiy 
-0.27 

1 63+ 017 
- L - D ' 3 -0.15 

1 64+ - 12 



00+" bb 
u - uu -o.oo 

77+0- 12 
u - ' ' -0.24 

70+ - 16 
'"-0.17 



ucll -0.31 
, ,^- + 0.20 
1 " !o -0.20 

1 09 + 015 
1 ' ua -0.15 



4 l 6 + l ti " 
*- iD -1.47 

4 30+ 45 
3 63+ 16 



ACIS-S 1' 
EPIC l' 
EPIC 4' 
RGS 1' 

RGS 3.5' 



1.39 
1.57 
1.39 



1.46 
1.65 
1.43 



„ 7 ,+0.002 
u - ' D1 -0.005 
„ -,,, + 0.002 
u - ,a ^-0.002 

771 + 0- 010 
' ' 1 -0.011 
„ 7 „„+0.007 

u. (a^_ 007 

2 + 0.035 

u - D1 ^-0 .036 

583+" 020 
u.ooa_ n 020 

„ -,,+0.063 
0.532_ Q89 

307+ 114 
U-JUl -0.032 



, , 7 +0.55 
-0.20 
, r, + 0.48 

n 8 ,+o.06 

0.»-i_o.06 
72+O.02 



o.»fr 



0.07 
+ 0.36 
0.21 

44+0-OU 
u, ^^-0.08 
0.37+0.04 



, +0.32 
-0.42 



112 
630 
931 
325 
552 



1.421 
1.408 
1.506 
1.335 
1.258 



0.48 
0.58 
1.75 
0.71 
1.66 



O80 + "- 11 80 + "- 10 44 + "- Is ()7I+ U1 , ■«+»■■« 
0.S0_ 10 0.80_ 09 0.44_ 18 0.71_ Q 15 2.24_ Q 36 



NGC5044 



EPIC l' 
EPIC 4' 
RGS l' 
RGS 3.5' 



0.02 
+ 0.09 
0.07 
+ 0.05 
0.04 

1 '°' ! -0.25 



0.65 
0.59 



0.42 
0.46 



0.04 
+ 0.10 
0.08 
+ 0.07 
0.06 



94+0- 05 
u - Mq -0.05 

n SQ+O- 33 

°' 89 -0.30 

49+ 14 
"' 4J -0.14 



89+ 04 
u.oy_ 04 



, + 0.19 
°-0Jj 



"-.6 + 0- 08 
O-OD_ ,08 



3 48+ 15 
J -*°-0.15 



477 1.118 

832 1.963 

295 1.262 

668 1.210 



1.83 
8.80 
4.60 
10.25 



IC1459 



EPIC l'no core 1 
EPIC 4' no corc-i 
RGS 1' 
RGS 3.5' 



0.39 
1.18 
4.96 



+ 0.14 
0.09 
+ 3.40 
0.84 
+ 0.04 
4.66 



, ,, + 2.68 
-1.79 
1 84+ 316 



171 0.828 

462 1.049 

66 1.105 

138 0.994 



0.90 
1.96 
2.02 
2.86 



Note: for RGS data, the Mg and Si 



"■Instruments used and their apertures. EPIC:joint fit for MOSl,MOS2,& pn data; RGS: joint fit for RGS1 & RGS2 data. EPIC source size: 
diameter of 1 and 4 ; RGS source size: including 90% and 99% of psf, around 1 and 3.5 in the cross-dispersion direction. One phase model: 
wabs(vapec+power) for EPIC, and wabs(vapec+power) *rgsxsrc for RGS. 

b One phase temperature in unit of keV. 

c Elcmcntal abundances of Fc, O, Mg, Si, Nc, S, and Ni relative to their Solar values (Grevess et al, 1! 
abundances were allowed to fit when thawing them can improve the fits, otherwise they were tied to Fe. 

d Degrees of freedom. 

°Rcduccd x 2 °f the fit. 

^Luminosity of the power-law component in unit of 10 4 "ergs/s and in the range of 0.33 keV to 2.48 keV. 
g Luminosity of the thermal component (vapec) in unit of 10 40 er<?s/s and in the range of 0.33 keV to 2.48 keV 
^Luminosity ratio of the thermal component to the power-law one. 

"We excluded the Fc XVII line at 15 A (i.e. excluding 14.75 — 15.40 Aj to correct the optical depth effect. 
J Wc excluded the central 13 radius region to eliminate the central non-hot-gas emission. 



Table 3. Two Phase Model Fitting for XMM EPIC and RGS spectra. 



Galaxy 



Camera/Size 



Tl(keV) T2(koV) 



Mg 



Power Vapecl Vapec2 Vapec/Power 



NGC720 



NGC1399 



NGC3923 



NGC4406 



EPIC l' 
EPIC 4' 



217+ 0161 
u - is -< -0.079 

n noo + 0-031 



0.567] 



+ 0.019 
-0.017 
562+°- ul " 



0.64 
1.10 



EPIC 1' 
EPIC 4' 
RGS 1' 
RGS 3.5' 



793+" :i » l:i i 71 <j+».0H5 
u - lad _(i nn/i , y — 0.070 



+ 0.38 
-0.17 
+0.11 
0-26 



0.30]; 

0.45 4 



EPIC l' 
EPIC 4' 



EPIC 1' 
EPIC 4' 
RGS l' 
RGS 3.5' 



EPIC 1' 
EPIC 4' 
RGS l' 
RGS 3.5' 



EPIC l' 
EPIC 4' 



EPIC l' 
EPIC 4' 
RGS l' 
RGS 3.5' 
RGS corrected l' 0.336 
RGS corrected 3.5 



-0.004 

n Kna+ - 002 

°- 809 -0 002 

n »nn+ 013 
0-800_o ooe 

660+ ' 052 

244+' U,4:1 

u - z ^_0.039 

247+° 024 
"""'-0.025 
-7-— FTFTTTTr 

U.OoU_ .023 
K-13+0 012 

017+ 026 
U.D-L I _o_034 

°- 2b6 -0.014 

7 , 4 +o.ooy 

u. '-3^_o.oo5 
„ 77 „ + 0.006 
u. I 'o_ .006 

0.682j:S;ol9 

656+° 126 
"■ DJD -0.055 

261 + "-™ 2 

U -^ D1 -0.139 

353+ - 029 
„ 33 6 + 0-013 



1 440+ 031 
1 '^ u -0.034 



0.555 



0.984±° 
0.938+° 
0.984 
0.702+° 
1.277+Jj 
1.313+.° 
1.021+° 
0.967+.° 



0.009 



TT7 



Trrr 

010 
1)09 
1109 
TTTV 

(KM 
1)20 
924 

023 
924 
TTT2- 

055 
923 
922 
-190 
29 4 
13S 
53 4 



1 70+0-Os 
L - '"-0.12 
. „,+0.05 
1 - JZ -0.05 



91 + 013 
u - J -0.09 
+0.1)7 



0.81 
0.72 



-0-11 
+ 0-25 
-0-19 



+0.24 
-0.12 
+ 0-22 
0.14 

"■ '"-0.09 

61 + °'° 5 
u - D1 -0.05 
, ,,+0.32 
^ -0.15 
A ^o + 0.12 
°- 63 -0.08 
-+0.12 



n t-,+0.46 
0.56_ 23 

76+ ' 73 
"• '°-0.32 
-, oo + O-ll 



0.49 
0.97 



+ 0.47 
0.33 
+ 0.93 
0.45 



. + 0.69 
D -0.70 



196 1.116 
500 1.094 



1.36 





+ 0.07 
0.07 



1.75"' 



-+T7779- 

-0.11 

1 48+°-° 5 
1 ' w -0.05 



, + 9.19 

-0.19 



l.ll 1 



+ 9.12 

0.13 
1 03+°° 6 



4 48+ 1 '- 34 
^■ qo -0.34 
+ 0.19 
' 0-10 



3-261 



0.27] 
0.24+ 



0.87 
0.33 
0.59 



1.55 
1.75 



. 09 
+ 0.02 
0.04 
+ 0.03 
0.05 
+ 0.11 
— 0.11 
+ 0.09 



-0.08 
■0.09 
-0.06 

u -^ u _ n in 



-0.10 
+ 0.05 
■0.05 
+ 0.10 
■0.08 

20+ - 04 
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NOTE: Table is the same as previous one, only that we use two phase model here: wabs(vapec-hvapec-hpower) for EPIC, and wabs(vapec+vapec+power) *rgsxsrc for RGS. For RGS 

second temperature will be frozen to the EPIC fits if it cannot be constrained. 
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